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Abstract 

We theoretically study multiphoton ionization of the Ca atom irradiated by the second (photon 
energy 3.1 eV) and third (photon energy 4.65 eV) harmonics of Ti:sapphire laser pulses (photon 
energy 1.55 eV). Because of the dense energy level structure the second and third harmonics of 
a Ti:sapphire laser are nearly single-photon resonant with the AsAp 1 P° and 4s5p 1 P° states, 
respectively. Although two-photon ionization takes place through the near-resonant intermediate 
states with the same symmetry in both cases, it turns out that there are significant differences 
between them. The photoelectron energy spectra exhibit the absence/presence of substructures. 
More interestingly, the photoelectron angular distributions clearly show that the main contribution 
to the ionization processes by the third harmonic arises from the far off- resonant 4s4p 1 P° state 
rather than the near-resonant 4s5p l P° state. These findings can be attributed to the fact that 
the dipole moment for the 4s 2 1 S e - As5p 1 P° transition is much smaller than that for the 4s 2 1 5 e 
- 4s4p 1 P° transition. 
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I. INTRODUCTION 

Above-threshold ionization (ATI) [1] is a process in which atoms absorb more than the 
minimum number of photons required to ionize and the photoelectron energy spectrum 
(PES) consists of a series of peaks that are equally separated by the photon energy. Now 
the ATI and multiphoton ionization (MPI) processes have been well-studied (2, 3|, especially 
for rare gas atoms. Recently, we have presented several new interesting features in the PES 
of a light alkaline-earth- metal atom, Mg, interacting with a short laser pulse 0, 0]. We 
have had a close look at the intermediate ATI peaks appearing in the PES and clarified 
their origin as the off-resonant excitation of the bound states. In this paper we extend our 
previous works on Mg to another alkaline-earth-metal atom which is less investigated in the 
literature, the Ca atom. 

During the last 30 years many theoretical and experimental investigations have been 
performed for Ca to obtain the atomic data and to understand its interaction with laser 
fields for the single-photon processes. The first extensive theoretical studies for the singlet 
series 1 S and l P of Ca were performed by Fischer and Hansen p] with a multi-configuration 
Hartree-Fock (MCFH) method which included correlations between the valence electrons. 

h 

Later on, Mitroy [7| calculated the energy level and oscillator strength (OS) for the low- 
lying levels of Ca using the frozen-core Hartree-Fock (FCHF) approach including a model 
potential. Aymar and co-workers |8j| extensively studied light alkaline-earth-metal atoms 
using the multichannel quantum defect theory and eigenchannel R-matrix combined with 
polarization potentials. Hansen and co-workers [9] used a configuration interaction (CI) 
approach, based on the B-spline basis functions and a model potential with dielectronic 
core polarization potential, to calculate term energies and wave functions for the singlet and 
triplet states of Ca. Recently, Fischer and Tachiev [lOj calculated the energy levels, transition 
probabilities, and lifetimes for the singlet and triplet spectra of Ca using the MCFH method 
with lowest-order relativistic effects included through a Breit-Pauli Hamiltonian. 

As for the MPI processes Benec'h and Bachau [HI] calculated the one-, two-, and three- 
photon ionization cross sections of Ca from the ground state with B-spline basis functions 
and CI procedure, where the FCHF approach and polarization potentials were employed to 
construct the atomic basis. Regarding the experimental MPI processes of Ca, there are only 
several works in the literature, all of which involve only ns and ps laser pulses. DiMauro and 
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co-workers |l2j investigated single and double ionization of Ca by 10 ns Nd:YAG laser pu 
at 532 and 1064 nm in the intensity range of 10 10 — 10 13 W/cm 2 . Shao and co-workers 
analyzed single and double ionization of Ca by 35 and 200 ps Nd:YAG laser pulses at 532 
and 1064 nm in the intensity range of 10 10 — 2 x 10 13 W/cm 2 . Lately, Cohen and co-workers 
experimentally and theoretically studied the two-photon ionization spectra of Ca in the 
374-323 nm wavelength range with a ~5 ns dye laser. Although theoretical data for MPI of 
alkaline-earth-metal atoms such as Ca by fs laser pulses would provide additional valuable 
information for the purpose of understanding the multiphoton ionization dynamics, such 
data are still missing in the literature. 

The aim of this work is to extend our previous investigations for MPI of Mg 0, [|], Q, [3] 
to the Ca atom: In this paper we study the MPI processes of Ca by the second (photon 
energy 3.1 eV) and third (photon energy 4.65 eV) harmonics of Ti:sapphire laser pulses. 
For this purpose we use a nonperturbative method to solve the time-dependent Schrodinger 
equation (TDSE) with two active electrons. In Sec. [IT] we construct an atomic basis set in 
terms of discretized states and use it in Sec. II III to numerically solve the TDSE. Numerical 
results such as the ionization yields, PES, and photoelectron angular distribution (PAD) are 
presented in in Sec. IIV1 Atomic units (a.u.) are used throughout this paper unless otherwise 
mentioned. 



II. ATOMIC BASIS STATES 



The Ca atom is a two-valence-electron atom with a closed ionic core Ca 2+ (the nucleus 
and the 18 inner-shell electrons ls 2 2s 2 2p 6 3s 2 3p 6 ) and the two valence electrons, 4s 2 l S e . 
Since it is a heavier alkaline-earth- metal element than Be and Mg and hence the core, Ca 2+ , 
is softer, there are more complexities as well as subtleties in the Ca atomic structure. As it 



is already mentioned in the literature 



13, 



181 ] there are several different approaches to solve 



the Schrodinger equation to describe the interaction of a one- and two- valence-electron atom 
with a laser field. Since the general computational procedure has already been presented 



in Refs. IsJ, uji Il9l-l2lll and the specific details about the atomic structure calculation of 

n n n 

Ca have been reported in recent works [9|, [Hi, 1221 ] . we only make a brief description of the 
method we employ. The field- free one-electron Hamiltonian of Ca + , H a (r), is expressed as: 



where V e ff(r) is the effective potential acting on the valence electron of Ca + , r represents 
the position vector of the valence electron, Z is the electric charge of the core (Z — 2 for 
a two- valence-electron atom), and I is the orbital quantum number. In our approach the 
effective potential, V e ff(r), consists of the FCHF potential (FCHFP) and the additional 
core-polarization potential which will be introduced in the next subsection. 

A. One-Electron Orbitals: Frozen-Core Hartree-Fock approach 

To describe the ionic core Ca 2+ we have introduced the effective potential in Eq. ([I]), 
which is a sum of the FCHFP, V l HF (r), and the core-polarization potential, Vfi^r, a s , n): 



V eff (r) = V l HF (r) + Vnr,a s ,r l ). 



(2) 



V t p (r, a s ,ri) describes the interaction between the ionic core and the valence electrons and 
can be written in the form of 



a. 



lf(r,a fl ,rj) = -— 1-exp 



2r 4 



-(r/n) 6 



(3) 



where a s is the static dipole polarizability of Ca 2+ and ri the cutoff radii for the different 



atained by 
M for the 



orbital angular momenta, I = 0,1,2,..., etc [19[. The values of r/ have been o 
performing the fittings of the one-electron energies to their experimental values 
four lowest states of s, p, d and / series of Ca + . We have used the following set of cutoff 
radii, tq = 1.5457, t\ = 1.5857, r2 = 1.8771, and r/>3 = 1.5530 together with the static 
dipole polarizability of Ca 2+ which is a s = 3.16 24|. The relatively large value of the cutoff 
radius for I — 2 is due to the fact the 3d orbital penetrates the ionic core much more than 
the other orbitals, and therefore the core-polarization potential is more sensitive for the 
I = 2 orbital than for other orbitals. We note that other theoretical papers 
slightly different values for the static dipole polarizability and cutoff radii. 
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With the aid of the mathematical properties of the B-spline polynomials 18J, Il9] to 
expand the one-electron orbitals, solving the one-electron Schrodinger equation for the non- 
relativistic one-electron Hamiltonian given by Eq. ([T]) is now equivalent to an eigenvalue 
problem. 
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B. Two-electron states 



The field-free two-electron Hamiltonian, H a (x\, r 2 ), can be expressed as 



(4) 



i=l 



where H a (ri) represents the one-electron Hamiltonian for the ith electron as shown in Eq. 
(PQ), and V(ri,r 2 ) is a two-electron interaction operator which includes the static Coulomb 
interaction l/|i"i — r 2 | and the effective dielectronic interaction potential [ij], 25]. Here 
i - ! and r 2 are the position vectors of the two valence electrons. To solve the two-electron 
Schrodinger equation for the Hamiltonian given in Eq. the two-electron states can be 
constructed within the CI approach. Namely, we use a linear combination of the products 
of the two one-electron orbitals to represent a two-electron states and diagonalize the two- 



electron Hamiltonian given in Eq. (j^J). This is so-called a CI procedure 
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421). 



For Ca, which is heavier than Be and Mg but still relative^ 

7-9 



y light alkaline-earth-metal 



ll| and hence it is sufficient 



atom, the LS coupling is known to give a fair description {7- 
to label a two-electron state by the following set of quantum numbers: principal, orbital, 
and spin quantum numbers for each electron, niUsi (i = 1,2), total orbital momentum L, 
total spin 5, total angular momentum J, and its projection M on the quantization axis. 
After the CI procedure the two-electron states may be most generally labeled by the state 
energy and the quantum numbers (L, 5, J, M), and furthermore the above state labeling can 
be simplified to (L,M) for singlet states (5 = 0). Having obtained the two-electron wave 
functions we can calculate the dipole matrix elements as well as OSs for both LP and CP 
laser pulses. 



III. TIME-DEPENDENT SCHRODINGER EQUATION 

By making use of the two-electron states which have been constructed in Sec. [TTJ we can 
solve the TDSE for the two-electron atom interacting with a laser pulse. The TDSE reads 

i JU(r 1; r 2 ; t) = [H a (r u r 2 ) + D(t)) *(r 1} r 2 ; t), (5) 

where ^/(r!,r 2 ;t) is the two-electron wave function for the two electrons located at I - ! and 
r 2 at time t, and H a (ri, r 2 ) is the field-free two-electron Hamiltonian shown in Eq. (jlj). The 
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time-dependent interaction operator, D(t), between the atom and the laser pulse is written 
in the velocity gauge and dipole approximation as, 

£>(*) = -A(t).( Pl +p 2 ), (6) 

where pi and p 2 are the momenta of the two electrons and A(t) is the vector potential of 
the laser field which is given by 

A(t) = A f(t) C os(ut). (7) 

In the above equation u is a photon energy and A = A e q is an amplitude of the vector 
potential with e q being the unit polarization vector of the laser pulse. The unit polarization 
vector is expressed in spherical coordinates, and q — 0, 1, and —1 correspond to the LP, 
right-circularly polarized, and left-circularly polarized fields, respectively. f(t) represents 
the temporal envelope of the laser field which is assumed to be a cosine-squared function, 
i.e., f(t) = cos 2 {irt/2r) where r is the pulse duration for the full width at half maximum 
(FWHM) of the vector potential A(£). The temporal integration range of TDSE in Eq. (J5j) 
is taken from — r to r. 

In order to solve Eq. (jSJ), the time-dependent two-electron wave function, \I/(ri, r 2 ; t), is 
expanded as a linear combination of two-electron states ^(ri, r 2 ; E n ): 

*(n,r 2 ;*) = ]T C EnLM (tWr h r 2 ;E n ), (8) 

n,L,M 

where C e„lm (t) is a time-dependent expansion coefficient for a two-electron state with an 
energy, E n , an angular momentum, L, and its projection on the quantization axis, M. Now, 
by substituting Eq. (jSJ) into Eq. (JSJ) we obtain a set of first-order differential equations for 
the time-dependent expansion coefficients Ce u lm (t) which reads 

i-T.C E n LM(f) = ^2 [En5 n n'5LL>5MM> — D n LAln>L>M'(t)]CE n/ L'M>(t), (9) 
dt n',L',M' 

where D nLMn i L > M /(t) represents the dipole matrix element between two singlet states defined 
by the quantum numbers (nLM) and {n'L'M'). This means that we have neglected the 
spin-forbidden transitions between the triplet and singlet states, which is a reasonably good 
assumption for a light atom such as Ca. Specifically in what follows, we assume that the Ca 
atom is initially in the ground state, 4s 2 1 S e (M = 0), i.e., 
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\C e„lm (t = -t)\ 2 = 5 nA 5 L o5 M Q- (10) 

Once we have obtained the time-dependent expansion coefficients Ce u lm by solving Eq. 
()9]), the ionization yield, Y, photoelectron energy spectrum, dP/dE e , and photoelectron 
angular distribution, dP/d8, can be calculated at the end of the pulse from the following 
relations: 



Y = 1 — £ | C EnLM {t = +r) 

n,L,M(E„<0) 



(E e ) = £ \C EeLM (t = +r)\ 2 , 



L,M 



and 



dP 

~d~e 



dP 

IK 



^{-ife^ph + 1 P h (cos9) C EeLM (t = +r) 

L,M 



(11) 



(12) 



(13) 



where E e represents the photoelectron energy, P\ 2 are the Legendre polynomials, I2 is the 
orbital momentum of the photoelectron, and 9 is the angle between the electric field and 
the photoelectron momentum vectors. 5(E e ) is the total phase shift which is the sum of 
the Coulomb and short-range scattering phase shifts. The total phase shift, S(E P ), can be 



extracted from the asymptotic behavior of the photoelectron wave function 
at large distances r — > 00: 
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,I26|, 



^ki 2 (r) -> \ l — sin 



irk 



kr + - \n{2kr) - 1 2 tt/2 + 5(E e 

K 



(14) 



where k = (2£' e ) 1 / 2 represents the momentum of the photoelectron. Since we employ the 
discretized technique to describe the wave functions in a rigid spherical box, the photoelec- 
tron wave function vanishes at the edge of the box (r = R). This means that the following 
relation always holds: 

kR + — ln(2A;i?) — 12^/2 + 5{E e ) = mn, where m is an integer, (15) 



which enables us to calculate the total phase shift, S(E e 
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IV. NUMERICAL RESULTS AND DISCUSSION 



In this section we present representative numerical results for multiphoton ionization of 
Ca from the ground state by the second and third harmonics of fs Ti:sapphire laser pulses. 
For the numerical calculation we have found out that a spherical box of radius R = 500 
a.u. and the total angular momentum up to L = 9 with 1800 states for each L gives a good 
convergence in terms of the ionization yield and PES. A number of 402 5-spline polynomials 
of order 9 with a sine-like knot grid is employed. Note that all the numerical results reported 
in this paper are calculated for the 20 fs (FWHM) cosine-squared pulse in the velocity gauge 
unless otherwise stated. 

Before solving the TDSE in Eq. (jSj), however, we must perform several checks regarding 
the accuracy of the atomic basis for the singlet states of Ca. The level structure of the 
singlet states of Ca is presented in Fig. [TJ The first ionization threshold lies at E ion = 6.11 
eV relative to the ground state. In Table I we compare our calculated energies with other 
theoretical results [3, E| and the experimental data for the first ionization threshold and the 
first few low- lying singlet states for L = S,P,D, and F. The experimental data are taken 
from the database of National Institute of Standards and Technology (NIST) [23] and the 
energies are taken with respect to the double ionization threshold, Ca 2+ . From Table I we 
notice that our results by the FCHF method provide energy values as accurate as other 
theoretical results for the ionization threshold and the first few low-lying states. 

The next step is to check the accuracy of the wave function in terms of the OSs. Table II 



presents comparisons of the OSs for a single 
theoretical p, and the experimental data 



photon absorption we have calculated with other 



23, 



28|-|3ll for L = S, P, D, and F in both length 
and velocity gauges for the following single-photon transitions: 4s 2 1 S' e — > 4s(4 — 6)p l P° 
and 3d4p 1 P°, 4s4p l P° -> 4s(5 - 7) a l S e and 4p 2 l S e , AsAp 1 P° ->■ 4s(4 - 6)d l D e and 
Ap 2 l D e , As3d 1 D e 4s(4 - 6)p 1 P° and 3dAp l P°, and finally 4s3d 1 D e -»■ 4s(4 - 6)/ l F° 
and 3d4p l F° . Clearly our results on OSs are in good agreement with other theoretical 
and experimental data. Relatively large discrepancies appear in the velocity gauge for 
the 4s3d x D e — > AsAp 1 P° and 4s3d 1 D e — > 3d4p 1 F° transitions. This may be due to 
the insufficient accuracy of the 3d orbitals. Fortunately the discrepancies appear for the 
transitions with very small values of the OSs, and therefore could hardly influence the 
outcome of the TDSE calculations we report in this work. 



8 



After we have checked the accuracy of the constructed atomic basis, we can now proceed 
to perform the time-dependent calculations by solving Eq. (Q under various intensities 
and photon energies for both LP and CP laser pulses. In the following calculations we 
employ 1800 two-electron states for each angular momentum up to L = 9 and carry out the 
numerical integration of TDSE [Eq. flU])] using a Runge-Kutta method. 



A. Ionization by the second harmonic of a Tirsapphire laser 

In this subsection we investigate two-photon ionization of Ca by the second harmonic 
of the Tksapphire laser at the photon energy u — 3.1 eV, which is schematically shown 
in Fig. |2j As already mentioned Ca has a relatively dense level structure and for photons 
in the visible range it is quite easy to be near-resonance with some bound states. Indeed, 
the detuning from the AsAp l P° state is only 0.17 eV, and considering the large value of 
the dipole matrix element for the transition 4s 2 1 S' e — > AsAp 1 P° (Table II), we expect a 
significant enhancement in the ionization signal. 

Figure [3]^a) shows the ionization yield as a function of peak intensity for LP (solid) and 
CP (dashed) laser pulses. The slope of these curves is about 1.85 up to the peak intensity 
of / = 5 x 10 11 W/cm 2 , after which the saturation takes place. The calculated slope is a 
little bit smaller than the prediction by the lowest order perturbation theory (LOPT) which 
gives a slope of 2 for two-photon ionization processes. Figure [3Jb) shows the ratio between 
the ionization yield by CP and LP laser pulses, Ycp/Y^p, as a function of peak intensity. 
In the low intensity regime (/ < 5 x 10 11 W/cm 2 ) the ratio slightly decreases with peak 
intensity and at peak intensity of 5 x 10 11 W/cm 2 it is about 1.25, which is smaller than 
the ratio predicted by the perturbation theory for the two-photon ionization cross sections 



of one- valence-electron atoms, crcp/ a LP = 1-4 [32M35] . 

In Fig. H]we plot the PES for the LP (solid) and CP (dashed) laser pulses at the peak 
intensity of / = 5 x 10 11 W/cm 2 . The ATI peaks exhibit some small structures on the left 
as well on the right wings which are equidistantly separated by the photon energy. In order 
to clarify the importance of the near-resonant bound state AsAp 1 P° we solved the TDSE 
with the AsAp 1 P° state artificially removed during the numerical integration. In Fig. |5^a) 
we compare the PES (solid) with that calculated without the state AsAp l P° (dashed) in the 
atomic basis. As a consequence the ionization signal is almost 4 orders of magnitude lower 
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and almost all substructures on the left and right wings of the ATI peaks disappeared when 
4sAp l P° is removed. In order to see the population dynamics of near-resonant AsAp 1 P°, 
we plot the population of 4s 2 1 S e (circles) and AsAp 1 P° (squares) states in Fig. Mb) as a 
function of time. Rabi oscillations take place between them. 

Finally, we show in Figs. E^a)-(B^d) the PADs at the photoelectron energies corresponding 
to the first four ATI peaks (see Fig. [4]) by the LP pulse at the peak intensity of / = 5 x 10 11 
W/cm 2 . Different ATI peaks exhibit different PADs, since different partial waves make 
different contributions with different total phase shifts. This is the reason why the PADs in 
Figs. E](a) and Etc) and also Figs. Mb) and E](d) resemble each other, since the accessible 
continua belong to the same parity. We now have a closer look at Fig. [6^a). The PAD 
has a typical profile for two-photon ionization from an initial S state with one secondary 
maximum at 9 = 90° and two minima at 9 = 54° and 126°, respectively. Similarly the PAD 
shown in Fig. |6^b) has a typical profile for three-photon ionization from an initial S state 
which exhibits two secondary maxima and three minima. Figures |6](c) and M,d) present the 
PADs at the third and fourth ATI peaks shown in Fig. HI 



B. Ionization by the third harmonic of a Thsapphire laser 

In this subsection we study two-photon ionization by the third harmonic of the Ti:sapphire 
laser at the photon energy of 4.65 eV, which is schematically shown in Fig. [7J Although 
the detuning from the As5p 1 P° bound state is only 0.1 eV which is even smaller than 
the case of the second harmonic in the previous subsection, we do not expect any impor- 
tant enhancement in the ionization process because the dipole moment for the transition 
4s 2 l S e — > As5p 1 P° is much smaller than that for the 4s 2 1 S' e — > AsAp 1 P° (see Table II). 
Actually the OSs for the 4s 2 1 S e — > As5p 1 P° transition is anomalously small. This is in 
contrast with the more regular decrease of OSs for the 3s 2 1 S e — > 3snp 1 P° (n = 3 — 6) 



transitions of Mg 



16] . 



First, we present the ionization yield in Fig. M, function of peak intensity for the 

LP (solid) and CP (dashed) laser pulses. The curves for the LP and CP laser pulses look 
almost the same with a slope of 1.98 at the peak intensities lower than I = 2 x 10 13 W/cm 2 , 
which agrees very well with the LOPT prediction. This implies that the near-resonance 
with As5p X P° makes very little contribution to the ionization yield. In Fig. E^b) the ratio 
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between the ionization yield for CP and LP laser pulses, Ycp/Y^p, is shown as a function 
of peak intensity. We can see that the ionization yield by the LP and CP laser pulses are 
almost equal below the saturation intensity. 

Figure M shows the PES by the LP (solid) and CP (dashed) laser pulses at the peak 
intensity of I = 10 13 W/cm 2 . Note the large difference of the peak intensities we have 
employed for Figs. H] and IU which, however, results in the similar amount of the ionization 
yields. Interestingly, the PES exhibits more substructures around each ATI peak, labeled 
as (a)-(d) in Fig. [9j where substructure (b) is hard to recognize due to the overlap with the 
ATI peak. The substructures are equidistantly separated by the photon energy of 4.65 eV 
and appear for both LP and CP pulses. These substructures in the PES resemble those we 
have seen for Mg at u = 4.65 eV [4J. The fact that the substructures appear for both LP 
and CP pulses implies that the certain bound states accessible in both laser polarization 
could be responsible for them. The procedure we have employed to identify the origin of 
the substructures is very similar to that we have used for Mg in our previous paper |4j. By 
inspection we expect that the bound states 4snp 1 P° (n=4,5, and 6) and 3dAp 1 P° could 
generate such substructures. Since we propagate the TDSE on the atomic basis, we can 
easily check this speculation by solving the TDSE by artificially removing the particular 
state under suspect, and comparing the PES with the original one obtained by the complete 
atomic basis. In Figs. [Tor a)- fTor d) we show the comparisons of the PES calculated under 
the same laser parameters with those for Fig. |9j We show the results obtained after the 
removal of a particular bound state (dashed lines), namely (a) AsAp l P°, (b) 4s5p l P°, (c) 
As6p 1 P°, and (d) 3dAp 1 P° upon solving the TDSE, in comparison with the PES for the 
complete calculation (solid lines) with the complete atomic basis. When the 4s4p l P° state 
is removed [Fig. [TOT a)]. the substructures on the left-side of each main peak are reduced 
or disappear, as highlighted by the circles. Similarly, by removing the 4s6p 1 P° and 3d4p 
1 P° states in Figs. fTUTc) and ITOTd). another small spikes labeled as (c) and (d) in Fig. [9] 
disappear. As for Fig. ITOT b) by artificially removing the near-resonant state 4s5p 1 P°, the 
height of the main ATI peaks is only slightly reduced since the state 4s5p 1 P° brings a small 
contribution in the ionization process (see Table II). These comparisons indicate that the 
physical origin of the substructures labeled as (a)-(d) in Fig. 19 is quite similar to that we 
have already found for the singlet as well triplet states of Mg |4j, |5j. Briefly, the off-resonant 
bound states such as 4s4p 1 P°, 4s6p 1 P°, and 3d4p 1 P° and the near-resonant bound state 
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4s5p 1 P° are the origin of the substructures. 

Finally, we show in Figs. [TTT a)- fTTT d) the PADs for the first four ATI peaks in Fig. [9] by 
the LP pulse at the peak intensity of / = 10 12 W/cm 2 . The peak intensity is chosen to be 
low to avoid any undesired complications at higher intensities. Again, different ATI peaks 
result in the different PADs. In order to examine the influence of the intermediate bound 
states with the 1 P° symmetry on the PAD we make a comparison of PADs by artificially 
removing the 4s4p 1 P° (dashed), 4s5p 1 P° (dot-dashed), 4s6p X P (dot-dot-dashed), and 
3d4p 1 P° (dot-dotted) states. The results are shown in Fig. [12j It turns out that the PAD 
is very sensitive to the removal of the off-resonant 4s4p 1 P° state (see Fig. [7]), while the 
change is very little when other bound states, including the near-resonant 4s5p 1 P° state, 
are removed. 



V. CONCLUSIONS 



We have theoretically studied multiphoton ionization of Ca by linearly and circularly 
polarized fs laser pulses at the photon energies of 3.1 eV and 4.65 eV in terms of the 
ionization yields, photoelectron energy spectra, and photoelectron angular distributions. 
At the photon energy of u> = 3.1 eV, the ionization process is strongly enhanced due to 
the presence of the near-resonant 4s4p 1 P° state which has a large dipole moment from the 
ground state, and the ionization yield is about four orders of magnitude larger compared with 
a case without a resonance. The photoelectron energy spectrum hardly shows substructures, 
because any possible substructures are buried between the strongly enhanced ATI peaks. In 
contrast, at the photon energy of u = 4.65 eV, the photoelectron energy spectrum exhibits 
many substructures due to the real excitation of the near-resonant 4s5p l P° state and some 
off-resonant bound states such as 4s4p 1 P°, 4s6p l P°, and 3d4p 1 P°, ... etc. Interestingly, 
the far off-resonant 4s4p 1 P° state still makes a very large contribution to the ionization 
processes in this case, which can be most clearly understood in terms of the photoelectron 
angular distribution. 
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Table I. Comparison of the energies for the first ionization threshold and the 
first few low-lying states of Ca. The energies are expressed in eV with 



respect to the energy of Ca 2 





Present 


Theory [7] 


Theory [9] 


NIST [23] 


^Ca+(4s) 


-11.87199 


-11.87112 


-11.87179 


-11.87172 




-18.011 


-17.950 


-17.988 


-17.98488 




-13.859 


-13.842 


-13.857 


-13.85406 


Ea s q s lge 


-12.973 


-12.920 


-12.950 


-12.93991 


E 4p 2 l S e 


-12.796 


-12.739 


-12.793 


-12.80404 


-£/4s4p i po 


-15.065 


-15.063 


-15.071 


-15.05236 


E^sSp ipo 


-13.427 


-13.429 


-13.436 


-13.43073 


E^gQp Ipo 


-12.798 


-12.805 


-12.821 


-12.81733 


E^dAp 1 P° 


-12.512 


-12.522 


-12.543 


-12.53782 


E/^sSd 1 D e 


-15.217 


-15.230 


-15.286 


-15.27586 


E/lsM 1 D e 


-13.321 


-13.322 


-13.377 


-13.36048 


E ip 2 l D e 


-12.901 


-12.914 


-12.939 


-12.93626 


E4, s 5d 1 D e 


-12.658 


-12.650 


-12.666 


-12.66359 


E^dAp 1 F° 


-12.925 


-12.934 


-12.964 


-12.95882 


E 4s 4f ipo 


-12.728 


-12.723 


-12.737 


-12.73494 


-^4s5/ IF" 


-12.428 


-12.424 


-12.430 


-12.42978 


E^sQf IF" 


-12.258 


-12.255 


-12.259 


-12.25899 



16 



Table II. Comparison of the absorption oscillator strengths (in atomic units and length/velocity 
gauge) between the few representative bound states with 1 S e , 1 P°, 1 D e , and 1 F° sym- 
metries. 



4s 2 1 S e -» 


4s4p 1 P° 


As5p 1 P° 


AsQp 1 P° 


3dAp 1 P° 


Present 


1.921/1.756 


4.115/3.563[-3] 


2.677/2. 151[-2] 


7.496/5.056[-2] 


Theory [7] 


1.82/1.781 


1.08/2.30[-3] 


3.72/3.14[-2] 


7.27/6.66[-2] 


Theory [9] 


1.745 


1.95[-3] 


3.62[-2] 


6.46[-2] 


NIST [23] 


1.75 




4.32[-2] 


7.01 [-2] 


Experimental [28] 


1.75 


9.0[-4] 


4.1[-2] 


6.6[-2] 


4s4p 1 P° -> 


4s5s ^ 


4s6s 1 S e 


Ap 2 l S e 


4s7s 1 S e 


Present 


0.112/0.091 


2.189/3.536[-2] 


0.110/0.129 


1. 648/1. 655[-2] 


Theory [7] 


0.118/0.114 


0.909/2.26[-2] 


0.120/0.148 




Theory [9] 


0.128 


7.50[-3] 


0.116 


1.01[-2] 


Experimental [29] 




9.0[-3] 


0.114 


1.33[-2] 


4s4p 1 P° 


AsAd 1 D e 


Ap 2 1 D e 


As5d 1 D e 


As6d 1 D e 


Present 


0.229/0.198 


0.573/0.604 


0.299/0.298 


5.768/5.814[-2] 


Theory [7] 


0.193/0.164 


0.550/0.531 


0.283/0.293 


5.28/5.38[-2] 


Theory [9] 


0.206 


0.548 


0.265 


4.527[-2] 


Experimental [29] 


0.207 


0.58 


0.28 


4.4[-2] 


4s3ci 1 P e -> 


AsAp l P° 


As5p 1 P° 


AsQp 1 P° 


3dAp 1 P° 


Present 


1.14[-3]/5.936[-2] 


5.730/7.876[-2] 


6.112/6. 756[-2] 


5.378/5.171[-2] 


Theory [7] 


8.759[-4]/6.332[-2] 


5.503/7.619[-2] 


5.533/6. 179[-2] 


5.26/5.73[-2] 


Theory [9j 


1.46[-3] 


5.85[-2] 


6.568[-2] 


5.166[-2] 


NIST [23] 




4.9[-2] 


7.5[-2] 


7.6[-2] 


Experimental [30, .31] 


1.02[-3] 


5.985[-2] 


5.66[-2] 


6.76[-2] 


4s3ci 1 P e 


3dAp 1 F° 


AsAf 1 F° 


4s5/ 1 F° 


4s6/ 1 F° 


Present 


6.877/3.397[-2] 


0.140/0.117 


5.657/5. 161[-2] 


3.053/2.842[-2] 


Theory [7] 


9.25/5.45[-2] 


0.129/0.111 


4.99/4.56[-2] 


2.23/2.27[-2] 


Theory [9] 


7.144[-2] 


0.131 


5.57[-2] 


3.038[-2] 


NIST [23] 


9.8[-2] 


9.39[-2] 


7.6[-2] 


3.2[-2] 


Experimental [30] 


9.8[-2] 


0.138 
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FIG. 1: Relevant energy levels of Ca. 
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FIG. 2: Relevant energy levels for ionization at the photon energy of 3.1 eV. The energy detuning 
from the 4s4p 1 P° bound state is 0.17 eV. 
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Peak Intensity (W/cm ) 



FIG. 3: (Color Online) (a) Ionization yield as a function of peak intensity by the linearly (solid) 
and circularly polarized (dashed) laser pulses at the photon energy of 3.1 eV. (b) Ratio of the 
ionization yield by the circularly polarized laser pulse, Yep, to that by the linearly polarized laser 
pulse, Y LP . 
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Photoelectron Energy (eV) 



FIG. 4: (Color Online) Photoelectron energy spectra by the linearly (solid) and circularly polarized 
(dashed) laser pulses at the photon energy of 3.1 eV and the peak intensity of 5xlO n W/cm 2 . 
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FIG. 5: (Color Online) (a) Comparison of the photoelectron energy spectra at the photon energy of 
3.1 eV. Solid line represents the result obtained from the complete atomic basis, while the dashed 
line represents the result from the atomic basis without the AsAp 1 P° state when solving the time- 
dependent Schrodinger equation, (b) Comparison of the populations of the 4s 2 1 S e (circles) and 
4s4p l P° (squares) states as a function of time. The laser pulse is linearly polarized and the peak 
intensity is 5xlO n W/cm 2 . 
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6(deg) 6(deg) 

FIG. 6: (Color Online) Photoelectron angular distributions at the photon energy of 3.1 eV, corre- 
sponding to the (a) first, (b) second, (c) third, and (d) fourth ATI peaks in Fig. [H The laser pulse 
is linearly polarized and the peak intensity is 5xlO n W/cm 2 . 
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FIG. 7: Relevant energy levels for ionization at the photon energy of 4.65 eV. The energy detuning 
from the 4s5p 1 P° bound state is 0.1 eV. 
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FIG. 8: (Color Online) (a) Ionization yield as a function of peak intensity by the linearly (solid) 
and circularly (dashed) polarized laser pulses at the photon energy of 4.65 eV. (b) Ratio of the 
ionization yield by the circularly polarized laser pulse, Yep, to that by the linearly polarized laser 
pulse, Y LP . 
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FIG. 9: (Color Online) Photoelectron energy spectra by the linearly (solid) and circularly polarized 
(dashed) laser pulses at the photon energy of 4.65 eV and the peak intensity of 10 13 W/cm 2 . (a)-(d) 
represent the substructures. 
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FIG. 10: (Color online) Comparison of the photoelectron spectra at the photon energy of 4.65 eV 
when the (a) 4s4p 1 P°, (b) 4s5p 1 P°, (c) As6p 1 P° , and (d) 3d4p 1 P° states are artificially removed 
from the atomic basis when solving the time-dependent Schrodinger equation. The laser pulse is 
linearly polarized and the peak intensity is 10 13 W/cm 2 . 



27 




e(deg) 9(deg) 

FIG. 11: (Color Online) Photoelectron angular distributions at the photon energy of 4.65 eV, 
corresponding to the (a) first, (b) second, (c) third, and (d) fourth ATI peaks in Fig. [9l The laser 
pulse is linearly polarized and the peak intensity is 10 12 W/cm 2 . 
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FIG. 12: (Color Online) Comparison of the photoelectron angular distributions by a linearly po- 
larized laser pulse at the photon energy of 4.65 eV when the 4s4p 1 P°, 4s5p 1 P°, 4s6p 1 P°, and 
3d4p 1 P° bound states of Ca are artificially removed from the atomic basis when solving the 
time-dependent Schrodinger equation. The peak intensity is 10 12 W/cm 2 . 
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